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We study particle generation by a strong electric field in lattice QCD. To avoid the sign problem 
of the Minkowskian electric field, we adopt the "isospin" electric charge. When a strong electric 
field is applied, the insulating vacuum is broken down and pairs of charged particles are produced 
by the Schwinger mechanism. The competition against the color confining force is also discussed. 
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Quarks interact with not only gluons but also pho- 
tons. Because the electromagnetic interaction is much 
weaker than the strong interaction, the electromagnetic 
interaction has been neglected in most lattice QCD simu- 
lations. However, the electromagnetic interaction can be 
essential in several situations [H. For example, high pre- 
cision simulations including dynamical QED effects are 
recently possible. Such simulations can reproduce elec- 
tromagnetic properties of hadrons. Also, strong back- 
ground electromagnetic fields are important for QCD 
phenomenology. Even if the electromagnetic coupling is 
smaller than the QCD coupling, strong electromagnetic 
fields can drastically affect hadron properties. A strong 
external magnetic field is one hot topic in lattice QCD 
i- 

Compared to the case of magnetic fields, the imple- 
mentation of electric fields is not simple. Two kinds of 
electric fields are possible: the Euclidean electric field 
and the Minkowskian electric field. (In the following dis- 
cussion, we choose the axial gauge Ai{x) = 0.) 

For the Euclidean electric field, the Euclidean gauge 
field A4^{x) is introduced as 



(1) 



in the Wick rotation. This gauge field A^ [x) respects the 
U(l) symmetry in the Euclidean space. As understood 
from the relation to the chemical potential /i, the Eu- 
clidean electric field cannot describe particle generation. 
The Euclidean electric field has been used for calculating 
electric polarizabilities in lattice simulations fs^] . The an- 
alytic continuation is needed to obtain the physical result 
in the Minkowski space. 

On the other hand, for the Minkowskian electric field, 
the vector potential Aq{x) is independent of the trans- 
formation 



i5o -f- qAo(x) 



-di + qA(^{x) 



(2) 



The U(l) gauge symmetry is lost in the Euclidean space 
and it must be fixed beforehand. This electric field is 
exactly the original one in the Minkowski space. The im- 
plementation of the Minkowskian electric field is difficult 
in lattice simulations. The vector potential breaks the 
anti-Hermitian property of the Dirac operator, which is 



the notorious sign problem. The Minkowskian electric 
field has been applied only in (partially) quenched lat- 
tice QCD 0. Moreover, a constant Minkowskian electric 
field is impossible in a periodic box because the vector 
potential is a real number, as explained later. 

The sign problem of the Minkowskian electric field is 
completely the same as that of a chemical potential. As 
seen in Eq. ([2]), the vector potential Aq{x) is regarded as 
a coordinate-dependent chemical potential. This means 
that the same strategy is applicable to the Minkowskian 
electric field and a chemical potential. The sign problem 
occurs in the physical electric charge q — diag(q„,(jd) = 
diag(2e/3, — e/3). To avoid the sign problem, we adopt a 
specific choice of the electric charges. 
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diag f^- 



(3) 



which we call the "isospin" electric charge. This pre- 
scription is inspired by an isospin chemical potential 
113 = which does not cause the sign problem 

[5| . The complex phases of the fermion determinants are 
canceled out between the w-quark and d-quark sectors. 
Thus, the fermion determinant is semi-positive and the 
sign problem does not occur. 

For numerical simulations, we consider the external 
Minkowskian electric field in a finite box with periodic 
boundary conditions. To satisfy the periodic boundary 
condition, we divide a finite spatial size L to two regions; 
L > z > L/2 and L/2 > z > 0. We set the vector 
potential as 



. J+^^o(.-f) (f>.>0) 



All the spatial components are zero, Ai{x) — 0. This 
configuration is depicted in Fig. [1] In this configuration, 
the electric field is applied in the z-direction. 



E[z) = -dMz) 



-E, (|>z>0) 
-Eo {L>z>^). 



The voltage difference is F = EoL/2. Since the vector 
potential Ao{x) is not the gauge field in the Euclidean 
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space, a coordinate-independent constant shift is not al- 
lowed. The overall constant must be fixed by hand. In 
this configuration, we have imposed the condition 

J d'^xAoix) = 0, (6) 

which is equivalent to the charge neutrality in the total 
volume. 
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FIG. 1: The configuration of the vector potential Ao{z) and 
the electric field E{z). The boundary condition is periodic. 



is meson condensation. Charged mesons can be formed 
when the voltage difference exceeds twice the lightest 
charged meson mass. The charge generation occurs only 
above this threshold. Because the lightest charged me- 
son is a pion, the threshold is twice the charged pion 
mass. This is similar to the charged pion condensation 
at a finite isospin density, where the charged pion is con- 
densed in /i3 > ttItt [!]■ The second possibility is the 
deconfinement. The electric field draws the charged par- 
ticles. This force is oriented in the opposite direction to 
the color confining force. When the electric field over- 
comes the confining force, the charged particles can be 
separated. The threshold of the voltage is the neutral 
pion mass. This effect is sensitive to the system vol- 
ume. When the quark mass is small, the spatial size of 
the quark- antiquark pair gradually increases. Thus, the 
charge density gradually appears in a finite volume even 
if the electric field is smaller than the confining force, 
while this effect is suppressed in a larger volume. Note 
that, quantitatively, these thresholds can be shifted from 
the original pion mass by the effect of the electric field. 



On the lattice, the vector potential is introduced as real 
link variables uo{x) = e^^-^ol^) u'^^{x) = e-93'4o(x)^ 
For example, the Wilson Dirac operator with the vector 
potential is 

D = 1-A*^[(1-7.)T,+ + (1+7.)I^.-] 

i 

74)^0 'T^4+ + (l + 74)wor4_] (7) 

with [T^+]x,y = U^,{x)5,j,+f,^y and [T^,J\x^y = 
Uj^{y)5x-fi^y This is the same implementation as an 
isospin chemical potential [HI, 0| . 

What we expect in the Minkowskian electric field is as 
follows. When a strong electric field is applied, quark- 
antiquark pairs are produced from the vacuum by the 
Schwinger mechanism [?'] . The quarks and the antiquarks 
fiow along the electric field. However, we cannot ob- 
serve such a nonequilibrium process, either the Schwinger 
mechanism or the flowing electric current. In lattice QCD 
simulations, the observable must be in an equilibrium 
state. The equilibrium state in the electric field is the 
final state after the charged particles flow and stop. In 
a finite box, the highest and lowest voltage regions exist 
somewhere. Positive charged particles stop at the lowest 
voltage region and negative charged particles stop at the 
highest voltage region. As a consequence, a nonuniform 
charge density distribution appears. When the voltage 
difference increases, the charge density grows. 

In the confinement phase, we need to take into account 
color confinement. Charged particles are created by the 
Schwinger mechanism, but they cannot fiow freely due to 
the confining force. There are two possibilities to sepa- 
rate the charged particles in the confinement phase. The 
schematic figure is shown in Fig. [2] The first possibility 
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FIG. 2: The meson condensation and the deconfinement. 

If we can apply a constant electric field in a periodic 
box, we can generate a permanent electric current. In the 
case of the Minkowskian electric field, we cannot gener- 
ate a permanent electric current. In periodic boundary 
conditions, a constant voltage gradient is impossible and 
only the configuration like Fig. [T] is possible. We cannot 
make a circuit with a steady electric current. (In the 
case of the Euclidean electric field or a magnetic field, a 
constant field is possible in a periodic box Q. This is 
because the U(l) link variable is a complex phase fac- 
tor which has 27r periodicity. Actually, a constant vector 
current was measured in the lattice simulation with a 
magnetic field flo'|.) 

We performed the two-flavor full QCD simulation with 
the plaquette gauge action and the Wilson fermion ac- 
tion. The vector potential is included both in the va- 
lence and dynamical fermions. The gauge coupling is 
/3 — 5.32441 and the hopping parameter \s n — 0.1665. 
The lattice spacing is a ~ 0.13 fm and the pion mass 
IS TTItt — 0.4 GeV [11|. The spatial latt ice volume is 
— (12a)'^ with periodic boundary conditions. We used 
two temporal lattice sizes; = 12 for the confinement 
phase and Nr = 4 for the deconfinement phase. 
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First, we demonstrate the competition between the 
electric field and the color confining force. We con- 
sider a "charged" heavy-quark potential. The heavy- 
quark potential is extracted from the Wilson loop. When 
the heavy quarks have electric charges, the link variable 
Ui_i{x) is multiplied to the SU(3) link variable U^{x) as 



Wc=iT\{{U^{x)u^{x)} 

loop 



W^5C/(3) n (8) 

loop 



For the charged heavy-quark potential in a constant elec- 
tric field, we calculate the rectangular (i? x T) Wilson 
loop in the z-t plane. From Eq. (U), the expectation 
value of the Wilson loop is 



(W^c(i?,r)) = (W^5C/(3)(i?,T))et 
and the charged heavy-quark potential is 



EoRT 



(9) 



aeV = 0.00 i — o- 
aeV = 0.48 1— ^ 
aeV = 0.96 i — b- 




FIG. 3: The charged heavy-quark potential Vc{R)- 



VciR) = Vsui3)iR)-7:EoR- 



(10) 



The sign of the second term depends on the electric 
charge and the region {L > z > L/2 or L/2 > z > 0). 
We here choose the repulsive case. In quenched QCD, 
the SU(3) heavy-quark potential Vsu(3) (R) is indepen- 
dent of the electric field, and it is the so-called Cornell 
potential. Thus, 



Vc{R) 



^Eq^ R+^+ const. (11) 



The electric field suppresses the linear confining poten- 
tial. In full QCD, Vsu{3){R) can be modified by the 
electric field through the dynamical quark. In Fig. [31 we 
show the numerical result of the charged heavy-quark po- 
tential. The data are fitted by Eq. (fTTj) . Although this 
is the full QCD simulation, the charged heavy-quark po- 
tential is consistent with the quenched QCD form ([TT]) . 
At aeV = 0.96, the linear confining potential disappears 
because the electric field is almost the same as the string 
tension, a'^eEo/2 = 0.08 ~ a'^a. 

Next, we calculated the charge density 



_ 1 d\nZ 1 / 1 
e oAo{x) " ^ 



dD ' 
aAo(x), 



(12) 



This charge density is equivalent to the isospin density 
[5[. In Fig.m we show the charge density distribution in 
the deconfinement phase. The charge density was fitted 
by a combination of two linear functions. When the volt- 
age is applied, the charge density becomes finite in the 
Aq ^ region. The positive (negative) charge density 
means the appearance of the positive (negative) charged 
particles, i.e., M-quarks and d-quarks (w-quarks and d- 
quarks). The charge density increases when the voltage 
increases. 

In Fig.O we show the voltage dependence of the charge 
density at 2 = 0. We calculated at the voltage aeV = 0, 
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FIG. 4: The charge density distribution 713(2). 



0.12, 0.24, 0.36, 0.48, 0.72, and 0.96. In the deconfine- 
ment phase, the charge density grows monotonically be- 
cause the charged particles fiows freely. In the confine- 
ment phase, the interpretation of the result is nontrivial. 
If the charge density were generated by the pion conden- 
sation, the charge density would be zero in eV < 2mT,. 
The pion mass at F = is arn^n- ~ 0.26 [Il[. The charge 
density is finite in aeV < 2am-„ ~ 0.52. The charge 
density in eV < 2m7r is generated not by the meson con- 
densation but by the deconfinement in a finite volume. 
The charge density is finite even at aeV = 0.24, which 
is slightly below m^r. This is not surprising because the 
pion mass will be shifted by the electric field. In a smaller 
voltage eV <^ m.,^, the charge density is not generated. 
This is consistent with the expectation that the QCD 
vacuum is an insulator at zero temperature. 

To know the asymptotic behavior, we need to take the 
large system limit L — > 00. We can consider two types 
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FIG. 5: The voltage dependence of the charge density m at 
z = 0. 



of the large system limit. In the L ^ oo limit with 
fixed V (i.e., E 0), the electric field cannot induce 
the deconfinement. The charge density appears only in 
eV > 2rn7r- On the other hand, in the L ^ oo limit with 
fixed E (i.e., V — > oo), the system energy diverges. An 
infinite number of charged mesons are produced from the 
vacuum. The deconfinement occurs only in eE/2 > a. 

We have simulated the Minkowskian electric field by 
introducing the isospin electric charge. Similarly, we can 
adopt two-color QCD, which can also avoid the sign prob- 
lem. In these frameworks, it is possible to study QCD 
in electric fields. We have analyzed the charge particle 
generation by the Schwinger mechanism. Although we 
cannot observe the nonequilibrium process itself, we can 
study some nonperturbative aspects of the particle gen- 
eration. The particle generation is characteristic in the 
Minkowskian electric field because the vacuum is stable 
in the Euclidean electric field. 
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